Energy conservation for dynamical black holes 
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An energy conservation law is described, expressing the increase in mass-energy of a general black 
hole in terms of the energy densities of the infalling matter and gravitational radiation. For a 
growing black hole, this first law of black-hole dynamics is equivalent to an equation of Ashtekar & 
Krishnan, but the new integral and differential forms are regular in the limit where the black hole 
ceases to grow. An effective gravitational-radiation energy tensor is obtained, providing measures 
of both ingoing and outgoing, transverse and longitudinal gravitational radiation on and near a 
black hole. Corresponding energy-tensor forms of the first law involve a preferred time vector which 
plays the role for dynamical black holes which the stationary Killing vector plays for stationary 
black holes. Identifying an energy flux, vanishing if and only if the horizon is null, allows a division 
into energy-supply and work terms, as in the first law of thermodynamics. The energy supply can 
be expressed in terms of area increase and a newly defined surface gravity, yielding a Gibbs-like 
equation, with a similar form to the so-called first law for stationary black holes. 

PACS numbers: 04.70.Bw, 04.30.Db, 04.70.Dy 



Introduction. Ashtekar & Krishnan 0, S recently 
derived an energy-balance equation for dynamical black 
holes, expressing the increase in mass-energy of a black 
hole in terms of the energy densities of the infalling mat- 
ter and gravitational radiation. This was an important 
step in ongoing efforts to develop a theory of black-hole 
dynamics, as opposed to the textbook theory of black 
holes 0, la, 1 which mostly concerns either stationary 
space-times or event horizons, which cannot be located 
by mortals. The Ashtekar-Krishnan equation describes 
how a black hole grows, but does not include the physi- 
cally important limit where it is starved of nourishment 
and ceases to grow. This Letter summarizes some key re- 
sults of an approach which can be applied to black holes 
in any state. In particular, an energy-conservation equa- 
tion is found in several different forms, including a gener- 
alization of the Ashtekar-Krishnan equation. Full details 
are presented in a longer article [8(. 

The key geometrical objects are trapping horizons 
0- EH , which are hypersurfaces H (in space-time) foli- 
ated by marginal surfaces. A marginal surface is a spatial 
surface on which one null expansion vanishes, where the 
null expansions 9± may be defined by 9±d 2 A = L±(d 2 A), 
where d 2 A denotes the area element and L± the Lie 
derivatives along future-pointing null normal vectors l±. 
This is a surface where outgoing or ingoing light rays are 
just trapped, neither converging nor diverging. For the 
author's local definition of black hole as a future (0_ < 
on H, for 9 + = 0) outer (L^9 + < 0) trapping horizon, 
it was shown that, assuming the dominant energy con- 
dition, the marginal surfaces have spherical topology (if 
compact) and that H is either spatial or null at each 
point, with its rate of area increase being respectively 
positive or zero 0. This is similar to Hawking's so-called 
second law for event horizons, but for a physically locat- 
able horizon. 



In this Letter, the new results apply to any trapping 
horizon, and can therefore be ap plie d to black holes, 
white holes, traversible wormholes [Tj and cosmological 
horizons. In comparison, the Ashtekar-Krishnan formal- 
ism applies only to spatial trapping horizons. The under- 
lying aim of this framework for black-hole dynamics is to 
provide a local, physical understanding of black holes, 
which can be used to guide and interpret more detailed 
studies, including numerical computations, and so help 
to make better contact with astronomical observations. 

Method. Einstein gravity is assumed, with space-time 
metric g. Labelling the marginal surfaces by a coordi- 
nate x, H is generated by a vector £ = d/dx, normal to 
the marginal surfaces. The coordinate freedom here is 
just x i — ► x(x) and choice of angular coordinates, under 
which all the key formulas will be invariant. To study 
derivatives off H, one may construct a dual-null foliation 
[l2L IT^|: two families of null hypersurfaces, forming the 
wavefronts of ingoing and outgoing radiation, intersect- 
ing in the marginal surfaces. Such a dual-null foliation 
always exists locally, and for spatial H, it is unique. For 
null H, it is not unique, leaving some subtleties to be 
resolved. Labelling the null hypersurfaces by dual-null 



coordinates x , one has normal 1-forms 



-dx d 



such that g 
„+- — „-i 



) = 0, a normalization function 



g^ 1 {n + ) n~) < 0, and l± = —g~ 1 {n T )/g + ~ arc 
the normal projections of the evolution vectors d/dx ± . 
One can fix g^ = — 1, where = denotes evaluation on 
H, but it will be retained here, since it cannot generally 
be so fixed away from H . 

The area of the transverse surfaces S is A — § s d 2 A 
and the area radius R = ^ A/An is convenient. As a 
measure of the active gravitational mass-energy enclosed 
by a surface, it is useful to take the Hawking energy Ei| 
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where c = 1 units are adapted to the dual-null method. 

First law. The first result is an energy conservation 
law which will here be called the first law of black-hole 
dynamics: 

/ Z-ipd 2 A+ I wd 2 A^ (2) 
ox Js Js dx 

where (ip,w) are defined below. The equation follows 
from the Einstein equation, specifically the null focus- 
ing and cross-focusing equations @ for (L±9±, L T 6±), 
employing the Gauss-Bonnet and Gauss divergence the- 
orems on S. The Ashtckar-Krishnan derivation has some 
similarities, but uses a different decomposition of the Ein- 
stein equation, based on the conventional 3+1 formal- 
ism for spatial H . Dual-null foliations are generally less 
familiar, but much better adapted to radiation, which 
propagates outwards in advanced time profile in 

retarded time x~ at given angle. 

The 1-form %j) is an energy flux (or the dual vector 
is an energy-momentum density) and the function w is 
an energy density, and they are each composed of matter 
and gravitational parts, ip — ipm+ipg, w = w m +w g . The 
matter parts are defined in terms of the matter energy 
tensor T as 

w m = — trace T/2 (3) 
V> m = T-VR + w m VR (4) 

where V is the covariant derivative and the trace is in 
the space normal to S. The expressions have the same 
form as those in spherical symmetry |15| . in which case 
the gravitational terms vanish. The coordinate forms are 

w rn = -g + ~T + - (5) 
(V>m)± = 9 + ~T ±± L T R. (6) 

As simple examples in spherical symmetry, for a mass- 
less Klein-Gordon field 0, T±± = (d±<p) 2 are the energy 
densities of radiation propagating in the directions, 
with T_| = 0, whereas one finds the electric energy den- 
sity w m = £ 2 /8ir for an electric field £ = q/R 2 7 where 
q is the enclosed charge [l^. Then one can interpret 
j>5 £ • i>md 2 A as a rate of energy supply and § s w m d 2 A as 
a rate of work. 

The gravitational energy flux is found to have compo- 
nents 

(^ g )±=9 + -\W±\?L T R/i2nG (7) 

where the traceless bilinear forms a± — -LL±h — 9±h 
are the null shear tensors, h denotes the metric of S and 
_L projection by h. The expressions for ip g have a simi- 
lar form to expressions for the energy flux of transverse 
gravitational radiation in several limits where it is well 
defined: at infinity in an asymptotically flat space-time 
|l6lll7j . for linearized gravitational radiation in the short- 
wave approximation [y] , in cylindrical symmetry |18| and 
in a quasi-spherical approximation |l9t I20L l2l| . Indeed, 



ip reduces to the Bondi flux at null infinity. Thus one 
may physically interpret tp g as energy flux of transverse 
gravitational radiation and <f s £ • ip d 2 A as the rate of en- 
ergy supply to the black hole by the infalling matter and 
gravitational radiation. Finally one has 

w B = ICIVSttG (8) 

where the 1-forms C(±) = ~-L((« T • V)rj '■ )/g^ are nor- 
mal fundamental forms of S, and one uses £ = Qja for H 
with 9± = 0. For reasons best seen in a spin-coefficient 
formulation, uu g can be interpreted as the ener gy d ensity 
of ingoing longitudinal gravitational radiation 22] . 

In summary, the first law J2J) expresses the rate of 
change of black-hole mass as a sum of two terms, in- 
terpreted respectively as rate of energy supply and rate 
of work. The terminology is analogous to that of the 
first law of thermodynamics, which expresses the change 
of internal energy (heat) as a sum of heat supply and 
work. Heat flux vanishes in thermal equilibrium, and a 
similar property holds for black holes: the energy flux 
£ • tp vanishes if H is null, whereas the work density w is 
generally non-zero for H of any causal nature. This is a 
physically important distinction because (as mentioned) 
H is spatial for a growing black hole, but becomes null in 
the limit where it ceases to grow and reaches equilibrium. 

Integral forms. When the distinction between energy- 
supply and work terms is unimportant, the first law may 
be written simply as 

dE r dR l2 , 

-d-x^f s ^ dA ^ 

where e = w + £ • ip/(dR/dx) is the combined energy 
density, which can be divided into that due to the matter 
and the gravitational field in the obvious way, e = e m +e s . 
A corresponding integral form of the first law is 

[E]^ f e^d 2 Adx (10) 
Jh ox 

where [E] denotes the change in E along the horizon, 
from one marginal surface to another. This expression 
uses the generator- volume element d 2 Adx, which is reg- 
ular as H becomes null. Ashtekar & Krishnan used 
the proper- volume element d 3 v — y/g xx d 2 Adx, where 
gxx = s(£j£)i which yields the proper- volume form 

[E] S f id 3 v (11) 
Jh 

where e — e{dR/ dx) / -Jg X x is the proper energy density. 
This form is technically regular in the limit where H 
becomes null, where g xx — > 0, but less useful than the 
differential or generator-volume forms, since the proper- 
volume element vanishes, d 3 v — > 0, with e being gener- 
ally finite. Conversely, the integrand e(dR/dx) of the 
generator- volume form vanishes in the null limit. 

Energy-tensor forms. The horizon induces a preferred 
time vector \ which is the curl of R in the normal space: 
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X is orthogonal to R and the transverse surfaces, x ' dR = 
0, J_x = 0, has normalization g(XiX) = —g~ 1 {dR,dR) 
and becomes null on a trapping horizon, g{x,x) — 0- 
This generalizes the Kodama vector in spherical symme- 
try |15l l23|. In particular, x reduces to the stationary 
Killing vector for Schwarzschild and Reissner-Nordstrom 
black holes, so it can be regarded as playing the tradi- 
tional role of a stationary Killing vector even for dynam- 
ical black holes. 

One can also introduce the vector r which is dual to 
the horizon-generating vector £ in the normal space: r 
is normal to H, g(£,r) = 0, _Lr = 0, has normalization 
g(r, t) = —<?(£, £) and is regular in the null limit, becom- 
ing null itself, r — > £. Then the matter energy density is 
found to have a remarkably simple form: 



dR 

1 dx 



T( X ,r). 



(12) 



In spherical symmetry, one actually has dE/dx — 
AT(x, t) for any normal vector x and its orthogonal dual 



r. Dividing the integrated flux as [E] = [E] 
follows that the integrated matter flux is 



[E] m = [ T( Xl T)d 2 Adx. 

J H 



[E] g , it 



(13) 



This is reminiscent of the usual definition of energy in a 
stationary space-time as a volume integral, with x replac- 
ing the stationary Killing vector. Ashtekar & Krishnan 
gave a similar form, which here would be 



[E] m = I T(x, 

JH 



f)d 3 v 



(14) 



where f = t / ^Jg xx is the unit normal vector to H. This 
works for spatial trapping horizons, but not in the null 
limit, where a unit normal vector does not exist. How- 
ever, it reveals that the Ashtekar-Krishnan permissible 
vector fields are gauge-fixed versions (with g H = —1, 
£ + = — £~) of the more manifestly invariant x/(dR/dx). 
The Killing-like vector x seems to have fundamental im- 
portance for dynamical black holes, defining a preferred 
flow of time outside the black hole. 

The integrated flux due to gravitational terms can 
be written in a similar way by introducing an effective 
gravitational-radiation energy tensor in the normal 
space, with components 



e±± = 

e ±T = 



|ct ± ||732ttG 
-|C(±)| 2 /8^V 



(15) 
(16) 



Since these components are non-negative, satisfies the 
dominant energy condition, implying that the gravita- 
tional radiation carries positive energy. The generally 
non-symmetric nature of is curious, but does give it the 
correct number of components for the energy-momentum 
density of ingoing and outgoing, transverse and longitu- 
dinal gravitational radiation. Note also that (cr±, C(±)) 



each have the correct number (two) of independent com- 
ponents for describing the respective radiation. The iden- 
tification and neat division of these modes is another suc- 
cess for the dual-null method. 

Then the differential form of the first law is 

<f(T( X ,T) + Q(x,r))d 2 A (17) 
dx J s 

and the generator-volume form is 



[E] Si / (T( X , t) + 6(x, r))d'A dx. (18) 

JH 

If the trapping horizon is spatial, one can also write 



[E] = I (T( X , 

JH 



f) + e( x ,T))d 3 v 



(19) 



which is closest to the Ashtekar-Krishnan form, having 
identified x an d ©■ These three forms perhaps best illus- 
trate the nature of the first law as an energy-conservation 
equation, expressing the increase in the mass-energy E of 
the black hole due to the energy densities of the infalling 
matter and gravitational radiation. 

Gibbs-like form. The so-called first law for station- 
ary black holes 0, 0, IS ^ s actually analogous to 
the Gibbs equation rather than the first law of ther- 
modynamics; both are assumed independently in non- 
equilibrium thermodynamics, formulated as a local field 
theory ^ IH HI To find such an equation for 

dynamical black holes, one needs a definition of surface 
gravity n which generalizes the standard definition for 
stationary black holes. In spherical symmetry, there is 
a definition which has several desired properties 0, |2^| , 
including a relation x~ (VA%) = ±kx with the same form 
as the standard one. A simple generalization, modifying 
an earlier definition 123, is 



k = g + ~ {2L-9+ + 6+0-) R/ '4 



(20) 



for 6+ = 0. This has the desired property that it van- 
ishes where L-6+ = 0, which is the previously defined 
condition for degenerate trapping horizons 9] . Then the 
rate of energy supply satisfies 



1 f f)A 

i-^^GitTrJ 2 ^ < 21) 



This is a remarkable generalization of the k5A/8ttG term 
in the so-called first law for stationary black holes, or 
the K,(dA/dx)/8nG term found previously in spherical 
symmetry 15]. Thus the desired Gibbs-like equation is 
found simply as 



dE _ 1 
~dx = 8ttGA J s 



9A J2 A 

K-—d z A + <f> w 
ox J s 



OR 

dx 



d 2 A^. 



(22) 



Since stationary black holes have a Hawking tempera- 
ture Kh/2nk and an inferred entropy S = Ak/AhG, the 
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above result suggests that the same is true locally for dy- 
namical black holes, by analogy with the original Clau- 
sius concept of entropy, as previously argued in spher- 
ical symmetry 30] . Defining a geometric entropy flux 
vector ip = 2iripk/fiK and a geometric entropy supply 
dSo/dx — § s £ • ipd 2 A, a geometric entropy conservation 
equation dS/dx = dS /dx can also be obtained 

When comparing k with the standard definition of sur- 
face gravity for stationary black holes, one faces the prob- 
lem that the dual-null foliation is generally not unique 
for null H, and therefore neither is \ nor K - It becomes 
unique in spherical symmetry, where n reduces to the 
standard surface gravity l/4Gm for Schwarzschild and, 
non-trivially, \/ m 2 — q 2 JG(m+\/m 2 — q 2 ) 2 for Reissner- 
Nordstrom black holes [15|. For Kerr black holes, a dual- 
null foliation giving the correct x (stationary Killing vec- 
tor) and k should exist, but is not known explicitly, de- 
spite a recent effort [3]J . It should also be noted that the 
definition l|20|l of surface gravity is inequivalent (for spa- 
tial horizons) to that of Ashtekar & Krishnan, and that 
their generalization of the so-called first law also has a 
different form to the Gibbs-like equation (1221) . Similar re- 
marks apply to the approach of Booth & Fairhurst ■ 

Conclusion. The most practical result here is probably 
the identification of the effective gravitational-radiation 
energy tensor 0. A new generation of gravitational- wave 
detectors offer an entirely new window on the universe, 
and a considerable community is attempting to produce 
predictions of waveforms from cataclysmic events such 
as black-hole mergers, mainly by numerical simulations. 
Classically, gravitational radiation is well defined only in 
weak-field regimes, which may be outside the numerical 
domain; the radiation-extraction problem. However, the 
new O is defined on a black-hole horizon H (as located by 
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